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Abstract 



The charged particle multiphcities of two- and three-jet events from the reaction e^e^ — > Z'' — > 
hadrons are measured for TP decays to Ught quark (uds) flavors. Using recent theoretical 
expressions to account for biases from event selection, results corresponding to unbiased gluon 
jets are extracted over a range of jet energies from about 11 to 30 GeV. We find consistency 
between these results and direct measurements of unbiased gluon jet multiplicity from T and 
TP decays. The unbiased gluon jet data including the direct measurements are compared to 
corresponding results for quark jets. We perform fits based on analytic expressions for particle 
multiplicity in jets to determine the ratio r = N^/N^ of multiplicities between gluon and quark 
jets as a function of energy. We also determine the ratio of slopes, r*^^) = (dNg/d7/)/(dNg/d|/), 
and of curvatures, r*^^) = (d^Ng/d|/^)/(d^Nq/d?/^), where y specifies the energy scale. At 30 GeV, 
we find r = 1.422 ± 0.051, r^^) = 1.761 ± 0.071 and r^^) = 1.98 ± 0.13, where the uncertainties 
are the statistical and systematic terms added in quadrature. These results are in general 
agreement with theoretical predictions. In addition, we use the measurements of the energy 
dependence of and to determine an effective value of the ratio of QCD color factors, 
Ca/Cf. Our result, Ca/Cf = 2.23 ± 0.14 (total), is consistent with the QCD value of 2.25. 
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1 Introduction 



The mean charged particle multiphcity of a gluon jet has often been measured in the anni- 
hilation of an electron and positron to hadrons, e"'"e~ hadrons. The usual method (see for 
example is to select three-jet quark-antiquark-gluon qqg final states for which the events 

and individual jets are defined using a jet algorithm such as the Durham p or Lucius jet 
finder. The particle multiplicity of a jet determined with this technique is found to depend on 
which algorithm is employed. Therefore, these jets and the associated qqg events are called 
"biased." 



In contrast, theoretical calculations usually define gluon jet multiplicity inclusively, by the 
particles in hemispheres of gluon-gluon [gg) systems in an overall color singlet. Quark jets are 
defined analogously as hemispheres of quark- ant iquark [qq) systems. The hemisphere definition 
of jets yields results which are independent of a jet finder. Therefore, these jets are called 
"unbiased." Unbiased gluon jet multiplicity has so far been measured only in T ||^, |^ and 
Z° [0-ljn| decays, corresponding to jet energies of about 5 and 40 GeV, respectively. In contrast. 



unbiased quark jet multiplicity is easy to measure and has been determined at many scales (for 
a recent review, see |T^). 



It is of interest to measure unbiased gluon jet multiplicity at other scales. Such measure- 
ments would allow a test of recent predictions |]13|, from Quantum Chromodynamics (QCD) 
for the scale dependence of the multiplicity ratio r between gluon and quark jets and for the 
related ratios of slopes r*^^^ and of curvatures r^^^: 

r = ^ , (1) 



^(,) ^ d^N,/d^y 



dNg/dy 

dN./dy ' ^ ' 

' ~ d2Ng/d2y ' 

where and are the mean particle multiplicities of gluon and quark jets, ?/ = In (Q/A), Q 
is the jet energy and A is the QCD scale parameter. Recently, a method to extract unbiased 
gluon jet multiplicity from biased e"'"e~ qqg events was proposed ||T5|, [1^, extending earlier 



formalism ||T^. This method provides an indirect means to determine the particle multiplicity 
of gluon jets, in a manner which corresponds to the hemisphere definition of jets. By combining 
measurements of found from this method with the unbiased measurements for N^, the ratios 
r, r^^-* and r'-^-* can be determined at a variety of scales and used to test the corresponding QCD 
results in a quantitative manner. 

The purpose of the present study is to test this method [^, |16[| to obtain unbiased mea- 



surements of Ng and to apply the results to determine the ratios, eqs. (|1])-(D. The data were 
collected with the OPAL detector operating at the LEP accelerator at CERN. We determine the 
ratios ([3)-(0) a function of scale and compare the results to recent QCD predictions 



Previously published measurements of the scale dependence of r |T8|, |T9| and r^^^ fT^, [T^ were 
based on biased jet samples, whereas experimental results for the ratio of curvatures r^^^ have 
not previously been presented. We also use the results on the scale dependence of unbiased 
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gluon and quark jet multiplicities to derive a measurement of the ratio of QCD color factors 
Ca/ Cf- 



2 Theoretical Framework 



Analytic expressions for the mean particle multiplicity of e"'"e three-jet events, valid to the next- 
to-leading order of perturbation theory (NLO, also called MLLA |^), were recently presented 
in 0]: 

^qqg = ^qq{L,k^,Lu) + ^^gg{k±,Lu) , (4) 

Nggg = Ngg {Lgg, k±^Lu) + ^ Ngg (^X.Le) ' 

These results, along with the others presented in this section, are based on massless quarks. 
The reason for the two different expressions, eqs. (^) and (^), is that there is an ambiguity in 
the definition of the gluon jet transverse momentum when the gluon radiation is hard. fc_L,L„ 
and k±^Le are the transverse momenta of the gluon with respect to the quark-antiquark system 
using the definition of either the Lund {k±^Lu) llll or Leningrad {k±^Le) groups, see 



and references therein for further discussion. Note that due to QCD coherence, the multiplicity 
of a gluon jet produced in e+e" annihilations depends on its transverse momentum and not its 
energy (see for example fl^])- ^qqg is the particle multiphcity of a three-jet event selected using 
a jet algorithm (see below). N^^ and N^^ are the multiplicities of two-jet qq and gg systems, 
given by twice Ng and N^, respectively. The scales L, Lgg, k±^Lu and k±^Le are defined as follows: 

L - in(^) , (6) 

Lm = In(^) , (7) 

k±,Lu = 2 1n(^) , (8) 

k±,Le = 2 1n(^) , (9) 

Sqq = Pq-Pq , (10) 
P±,LU = , (11) 



P±,Le = , (12) 

V 

with s = -E^.m. ) ^ig — Pi 'Pg = 9) ^) 5 Pqi Pq Pg 4-momenta of the g, q and g. L specifies 
the e"'"e~ cm. energy (Ecm.) and Lgg the energy of the qq system in the qq rest frame. 

The multiplicity of the gg system in this formalism, N^^, depends only on a single scale: 
k±,Lu in eq. (^ or k±^Le in eq. (^. This dependence on a single scale is a statement that N^^ is 
unbiased, i.e. Ngg {k±) in eq. (§) or is equivalent to the inclusive multiplicity of a gg event 
from a color singlet source produced at the same scale k±, to NLO accuracy. In contrast, the 
multiplicity of the qq system, Ng^, depends on two scales: L and k±, 

,Lu in eq. (|^) or Lgg and 
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k±,Lu in eq. (j^). This dependence on two scales implies that the multiplicity of the qq system is 
biased, depending on a hard scale L or Lgq for the production of the quark and antiquark jets 
and on a cutoff k± ^ fc±,LM below which the gluon jet is not resolved. This accounting for the 
bias in quark jet multiplicity due to the jet finder criteria used to select the qqg events is the 
difference between the theoretical results used here WB and those presented previously 



The above expressions are valid for jet algorithms which employ a transverse momentum 
cutoff k± to resolve "two-jet" qq from "three-jet" qqg events. Examples of such algorithms are 
the Durham P], Cambridge |^ and Lucius jet finders. The formalism is not expected to 



be valid for the Jade or cone jet finders which use invariant masses or angles rather 



than transverse momenta to separate the two- and three-jet event classes. 



According to the prescription in fl^l, the resolution scale k± should be adjusted separately 
for each event so exactly three jets are reconstructed. In contrast, a fixed value of the resolution 
scale results in truncation of higher order radiation, biasing the gluon jet's properties. Note 
that any radiation ("sub-jet") emitted within the quark or gluon jet must necessarily have a 
smaller transverse momentum than the gluon jet itself: otherwise the role of the "gluon jet" 
and "sub-jet" would be reversed. Thus the transverse momentum of the gluon jet, k±^giuon, 
defines an effective cutoff for sub-jet radiation, i.e. k± = kj_^giuon ~ ^x.Lm- 



In the formalism of ||l5l, the biased quark jet multiplicities N^g (J, k±^Lu) in eqs. (§) and (|^) 
(with J representing either L or Lqg) can be obtained from the corresponding unbiased multi- 
plicities 'Nqq (J'), which depend only on a single scale J', through the relation: 

dN -( D 

N,, (J, k^,Lu) = N,, (JO + ( J - J') ' (13) 

where J' = k±^Lu + Cq with Cg = 3/2. The unbiased term N^^ (i^) is equivalent to the mean 
multiplicity of inclusive e"'"e~ hadrons events as a function of = Ecm.- 

Finally, a relationship between the scale evolution of unbiased gluon and quark jet multi- 
plicities is given ||T5[| : 

d (L') _CaL_ aoCr_\ dN,g(L) ^^^^ 



dL Cf \ L J dL 

where L' = L + Cg — Cq, Cg = 11/6, ao = 6 Ca/(11 Ca — 2 n/), Cr = 10 7r^/27 — 3/2, and n/ equals 
the number of active quark flavors, taken to be = 5. Ca and Cp are the QCD color factors, 
3 and 4/3, respectively. 



3 Detector and data sample 



The OPAL detector is described in detail elsewhere [p4t p5| . The tracking system consists 



of a silicon microvertex detector, an inner vertex chamber, a large volume jet chamber, and 
specialized chambers at the outer radius of the jet chamber which improve the measurements 
in the z-direction.Q The tracking system covers the region | cos6'| <0.98 and is enclosed by a 

^Our coordinate system is defined so that z is the coordinate parallel to the e^ beam axis, r is the coordinate 
normal to the beam axis, (fi is the azimuthal angle around the beam axis and 9 is the polar angle with respect to z. 
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solenoidal magnet coil with an axial field of 0.435 T. Electromagnetic energy is measured by a 
lead-glass calorimeter located outside the magnet coil, which also covers | cos 6*1 <0.98. 

The present analysis is based on a sample of about 2 283 000 hadronic annihilation events, 
corresponding to the OPAL sample collected within 0.3 GeV of the T? peak that includes 



readout of both the r-0 and z coordinates of the silicon strip microvertex detector ||25|. The 



procedures for identifying hadronic annihilation events are described in [2f: 



We employ charged tracks measured in the tracking chambers and clusters measured in the 
electromagnetic calorimeter. Charged tracks are required to have at least 20 measured points 
(of 159 possible) in the jet chamber, or at least 50% of the points expected based on the track's 
polar angle, whichever is smaller. In addition, the tracks are required to have a momentum 
in the direction perpendicular to the beam axis greater than 0.05 GeV/c, to lie in the region 
I cos 6'| < 0.96, to point to the origin to within 5 cm in the r-0 plane and 30 cm in the z direction, 
and to yield a reasonable per degree-of-freedom (d.o.f.) for the track fit in the r-0 plane. 
Clusters are required to have an energy greater than 0.10 GeV if they are in the barrel section of 
the detector (| cos 6*1 < 0.82) or 0.25 GeV if they are in the endcap section (0.82 < | cos 6*1 < 0.98). 
A matching algorithm is employed to reduce double counting of energy in cases where charged 
tracks point towards electromagnetic clusters. Specifically, the expected calorimeter energy of 
the associated tracks is subtracted from the cluster energy. If the energy of a cluster is smaller 
than that expected for the associated tracks, the cluster is not used. Each accepted track and 
cluster is considered to be a particle. Tracks are assigned the pion mass. Clusters are assigned 
zero mass since they originate mostly from photons. 

To eliminate residual background and events in which a significant number of particles is 
lost near the beam direction, the number of accepted charged tracks in an event is required to 



be at least five and the thrust axis |2^ of the event, calculated using the particles, is required 
to satisfy | cos(6'thrust)| < 0.90, where 6'thrust is the angle between the thrust and beam axes. The 
residual background to the sample of hadronic events from all sources is estimated to be less 
than 1% [H. 



4 Selection of uds events 



The theoretical formalism we employ is based on massless quarks, as stated in Sect. |^. There- 
fore, beginning with the data sample described in Sect. ^, we select light quark (u, d and s) 
events for our study. 

The uds event tagging is based on the signed impact parameter values of charged tracks with 
respect to the primary event vertex, dgign. The magnitude of dsign is given by the distance of 
closest approach of the track to the event vertex in the r-0 plane. The sign of d^ign is determined 
as follows. Jets are reconstructed in each event using the cone jet finder [^.0 For each track 
assigned to a jet by the jet finder, hemispheres are defined by the plane perpendicular to the 
corresponding jet axis. The sign of d^ign is determined using the crossing point in the r-0 plane 
between the track and jet axis to which it is assigned, dsign is positive if the crossing point is 

^The jet resolution parameters used are: a cone half angle of 0.55 radians and a minimum jet energy of 
5 GeV. 
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in the forward hemisphere and negative otherwise. The forward hemisphere is the hemisphere 
which contains the jet. Because of the relatively long lifetime of hadrons containing c and b 
quarks, the distribution of dgign is strongly skewed toward positive values for c and b events 
but not for uds events. 

Note that the jet finding procedure described in the previous paragraph is employed solely 
for the purpose of identifying uds events and not for the analysis of gluon jet multiplicity 
presented below in Sects. 

Charged tracks are selected for the uds tagging procedure if they have r-0 coordinate in- 
formation from at least one of the two silicon detector layers, a momentum of 1 GeV/c or 
larger, and a maximum distance of closest approach to the primary event vertex in the r-0 
plane of 0.3 cm with a maximum uncertainty on this quantity of 0.1 cm. In addition they must 
be assigned to a jet using the cone jet finding procedure mentioned above. If no track in an 
event satisfies these requirements (1.5% of the events), the event is eliminated. The number of 
tracks which meet these requirements and have dsignl crd^^gn >3.0 is determined, where o^d^^^^ is 
the uncertainty associated with dsign- An event is tagged as a uds event if this number is zero. 
In total 1 109 017 events are tagged. The number of events in our final event sample and the 
corresponding uds purity are presented in Sect. ^ 



5 Analysis method 

Three-jet events are defined using a jet algorithm. For our standard analysis we employ the 
Durham jet finder p. As a systematic check, we use the Cambridge |21[] and Lucius jet 
finders (see Sect. |^). We choose these algorithms because they employ a cutoff in transverse 
momentum to separate two- and three-jet events, in correspondence to the theoretical formalism 
of |]15|, In particular, this means that eqs. (H), (||), ([T3|) and (|I^ are expected to be 

applicable for these jet algorithms. 

The resolution scale of the Durham jet finder, ijcuu is adjusted separately for each tagged 
uds event so exactly three jets are reconstructed. For the Cambridge jet finder, the resolution 
scale is again Ucut- For Lucius, the corresponding parameter is djom- We apply this procedure 
rather than use a fixed resolution scale to avoid introducing a bias in the gluon jet properties as 
discussed in Sect. |^. The jets are assigned energies based on the angles between them, assuming 
massless kinematics (see for example p8[). The jets are ordered from 1 to 3 such that jet 1 has 
the highest energy. Events are retained if the angles between the highest energy jet and the 
other two are the same to within 3°, the so-called "Y events. A schematic diagram of a Y 
event is shown in Fig. |I|. For Y events, the scales L^^, k±^Lu and k^^Le (eqs. (0)-(i|)) depend only 
on Ec.m. and one inter-jet angle, conveniently chosen to be 9i (see Fig. |I|). For 35° <9i< 120°, 
the range of 6i we employ for our gluon jet analysis (Sect. |OD , 22 365 events are selected: this 
is our final event sample. We require 9i > 35° in order that the two lower energy jets be clearly 
separated from each other. Note that since 6*2 and 63 (see Fig. |l|) are not exactly equal, the 
energies of jets 2 and 3 are not exactly the same. 

The uds event purity of this sample is estimated using the Jetset Monte Carlo multihadronic 



^ Y events were first studied in pq, E9l ; the overall charged multiplicity of Y events was first presented in pi 
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Figure 1: Schematic representation of a three-jet quark- ant iquark-gluon qqg event with a Y 
event topology in which the angle between the highest energy jet and each of the 

two lower energy jets is about the same. The jets are ordered in energy such that jet 1 has 
the highest energy. The angle 9i opposite the highest energy jet is used to specify the event 
topology. In our analysis, the gluon jet is assumed to be jet 3. 



event generator including initial state photon radiation, detector simulation and the 



same analysis procedures as are applied to the data. We use a sample of about 6 000 000 Monte 
Carlo events generated with version 7.4 of the program and the parameters given in The 
estimated uds purity is found to be (78.5 ± 0.2 (stat.))%. The Monte Carlo predicts that 
70% of the background events are c events and 30% are b events. The ratio of the number 
of events in the final sample to the number in the initial multihadronic Z° decay sample is 
(0.980± 0.007 (stat.)) x lO'^ for the data and (1.036 ±0.004 (stat.)) x lO'^ for the Monte Carlo. 

For simplicity, we identify the gluon jet by assuming it is the lowest energy jet in an event, 
i.e. jet 3. For the majority of Y events, 6i << 62, 63- For these events, jet 1 has an overwhelming 
probability (typically about 97%) to be a quark (or antiquark) jet, while jets 2 and 3 have about 
an equal probability to arise from a quark or gluon. Since 62 (see Fig. |T|), the scales Lgg, 
k±,Lu and k±^Le are almost unchanged if the gluon is actually jet 2 rather than jet 3. Similarly, 
as the energy of jet 3 increases so that 61^62^0^, all three jets have about the same probability 
to be the gluon jet and the determination of scales is again insensitive to whether the gluon 
jet assignment is correct or not. This relative insensitivity to gluon jet misidentification is our 
motivation for employing Y events. 

The scales Lgg, k±^Lu and k±^Le and the three-jet event charged particle multiplicity N^^^ are 
measured as a function of 9i. Lgg, k±^Lu and k±^Le are determined using all accepted charged 
and neutral particles. For the QCD scale parameter A, used in the definition of these last three 
variables (also for the definition of L, see eq. (||)), we choose a value of 0.20 GeV. We verified 
our results are not sensitive to this choice. 



The measured distributions of N^^^, Lgg, k±^Lu and fc_L,Le are corrected for experimental 
acceptance and resolution, the effects of initial state photon radiation, and misidentification of 
uds events, using bin-by-bin multiplicative factors derived from Jetset. The corrections for Lgg, 
k±,Lu and k±^Le also account for gluon jet misidentification, i.e. for the contributions of events 
in which the gluon jet is not jet 3. The data are corrected to the hadron level. The hadron 
level does not include detector response or initial state radiation and treats all charged and 
neutral particles with lifetimes greater than 3 x 10~^° s as stable: hence charged particles from 
the decays of Kg and weakly decaying hyperons are included in the corrected distributions. 
The correction factors are determined as described in [Q. The corrections are typically 1.08 
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9i (degrees) 


AT eft. 


(GeV) 


P±,Lu (GeV) 


P±,Le (GeV) 


0-5 


9.31 ±0.13 


67.27 ±0 


18 


1.227 ±0.007 


1.595 ±0.016 


5-10 


14.739 ± 0.026 


66.182 ±0 


023 


2.104 ±0.003 


2.897 ±0.005 


10-15 


17.279 ± 0.024 


64.385 ± 


024 


3.291 ±0.003 


4.658 ±0.005 


15-20 


19.282 ±0.038 


62.679 ±0 


048 


4.885 ±0.005 


7.102 ±0.012 


20-25 


20.540 ± 0.056 


61.530 ±0 


070 


6.469 ± 0.009 


9.581 ±0.023 


25-30 


21.533 ±0.073 


60.945 ± 


086 


8.008 ±0.013 


11.975 ±0.035 


30-35 


22.417 ±0.092 


60.597 ±0 


097 


9.518 ±0.017 


14.315 ±0.048 


35-40 


23.23 ±0.11 


60.02 ±0 


11 


11.057 ±0.022 


16.790 ± 0.064 


40-45 


23.75 ±0.13 


59.61 ±0 


12 


12.553 ±0.027 


19.194 ±0.079 


45-50 


24.12 ±0.15 


59.04 ±0 


13 


14.090 ± 0.033 


21.751 ±0.098 


50-55 


25.24 ±0.17 


58.88 ±0 


14 


15.514 ±0.038 


24.01 ±0.11 


55-60 


25.49 ±0.19 


58.45 ± 


14 


16.991 ±0.043 


26.49 ±0.13 


60-65 


25.66 ±0.21 


57.94 ± 


15 


18.423 ± 0.049 


28.98 ±0.15 


65-70 


26.34 ±0.23 


57.32 ±0 


15 


19.922 ±0.055 


31.68 ±0.17 


70-75 


26.47 ±0.25 


56.55 ±0 


16 


21.443 ±0.061 


34.56 ±0.19 


75-80 


26.79 ±0.26 


55.81 ±0 


16 


22.879 ±0.065 


37.36 ±0.21 


80-85 


27.10 ±0.28 


55.89 ±0 


14 


24.018 ±0.064 


39.21 ±0.20 


85-90 


27.29 ±0.31 


55.04 ± 


14 


25.388 ±0.065 


42.04 ±0.21 


90-95 


27.37 ±0.31 


54.29 ± 


13 


26.687 ±0.066 


44.80 ± 0.22 


yo- iuu 


Z ( . ( Z IE u.oo 


53.48 ± 


12 


97 Q97 -(- n rifip; 


/I7 KQ _|_ n 99 


100-105 


28.42 ±0.36 


53.384 ±0 


097 


28.703 ± 0.054 


49.00 ±0.18 


105-110 


28.31 ±0.35 


52.779 ±0 


077 


29.620 ± 0.046 


51.15 ±0.15 


110-115 


29.08 ±0.36 


52.416 ±0 


053 


30.240 ± 0.032 


52.58 ±0.11 


115-120 


28.62 ±0.37 


52.328 ±0 


032 


30.497 ±0.020 


53.115 ±0.066 



Table 1: Measurements of the mean charged particle multiphcity N^^^ of three-jet uds flavor Y 
events from Z° decays, as a function of the angle 6i between the two lowest energy jets. The 
jets are defined using the Durham jet-finder. The results for the quark jet scale ^/s^ and the 
gluon jet scales p±^lu and p±^Le are also given. The uncertainties are statistical only. These 
data are shown in Fig. 0. 



for Ng^g, 0.98 for Lgg, 1.03 for k±^Lu, and 1.05 for k±^Le, and do not vary much with 9i. There is 
good agreement between the distributions of the uncorrected data and the Monte Carlo sample 
when the latter includes initial-state radiation and detector simulation, with the same analysis 
procedures applied to the simulated events as to the data. 

Corrected distributions of N^^^, y/s^, Pl,Lu and p^^Le are shown in Fig. ^ The last three 
variables are obtained from the corrected distributions of L^^, and k±^Le using eqs. (|^)- 

(|^). We prefer to present y^s^, p±,lu and p±,Le rather than L^,, k^^^iu and fc±,Le because they 
provide a clearer illustration of the physical scales in our analysis. These data are listed in 
Table ||. Besides the results from our standard analysis, the results found using the Cambridge 
and Lucius jet finders are shown in Fig. 0. For about 2% of the events, the resolution parameter 
Hcut of the Cambridge algorithm could not be adjusted to yield three reconstructed jets, see [|4 



for a discussion of this feature of the algorithm. These events were rejected for the Cambridge 
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Figure 2: (a) The mean charged particle multiphcity of three-jet uds flavor Y events from Z° 
decays, selected using the Durham 0, Cambridge and Lucius P] jet finders, as a function of 
the opening angle 9i; (b) the corresponding scales ^ys^, P±,lu and p±,Le, see eqs. (|TD|)-(|T2p. The 
data have been corrected for initial state photon radiation, detector response, misidentification 
of uds events, and (for part (b)) misidentification of gluon jets. The uncertainties (too small 
to be visible) are statistical only. 
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Ae//. (GeV) Normalization K (d.o.f.) 


Ng^- (L), uds quark jets 
{L), gluon jets 


0.190 ±0.032 0.136 ±0.009 6.1 (16) 
0.600 ±0.059 0.201 ±0.011 28.1(17) 



Table 2: Results of fits of the 3NL0 expressions for quark and gluon jet multiplicities to 
measurements of charged particle multiplicity in unbiased qq and gg events, see the solid curves 
in Figs. and The fits are performed assuming five active quark flavors. The uncertainties 
are statistical only. The P^r degree-of-freedom (d.o.f.) of the fits are also given. The 
values are determined using the statistical uncertainties of the data and not the systematic 
terms. 



jet finder based study. 

We invert eqs. (0) and (|) to obtain expressions for the unbiased gluon jet charged particle 
multiplicities N^^- (A;_|_): 

K (f'^^Luiei)) = 2\wi-(L,e^)-KUL,k^,Lu{ei))] , (is) 



99 
Tch. 



qqg V^) " i/ ^ ' qq 
Tch. I T Q \ Mc/i 



N^^- {k^,Le m = 2 [N-;^ (L, e,) - W^- (L,, (6,), k^,Lu (^i))J , (16) 
where N^^^, Lgg, k±^Lu and k±^Le are determined from the Y events as described above. 



To find the N^^- (L, k±^iu) terms in eq. (|T5D, we employ two methods. First, for the standard 



analysis, we perform a direct measurement. Specifically we determine the particle multiplicity of 
two-jet uds events from Z° decays (i.e. at a fixed value of L corresponding to -Ec.m. = 91.2 GeV) 
as a function of the jet resolution scale k±^Lu- This result is presented in Sect. ^A\ . Second, as a 
systematic check (see Sect. we evaluate the analytic expression, eq. ([T3|), as explained in the 
following paragraph. To find the N^^- {Lgg, k± Lu) terms in eq. ([iBD , we utilize only the second of 
these methods, i.e. the one based on eq. (0). A direct measurement is not straightforward in 
this case because Lgg, unlike L, is a variable quantity. A direct measurement of N^^' {Lgg, k±^Lu) 
at the scales relevant for our analysis would in principle require a determination of the particle 
multiplicity in two-jet events as a function of p±^Lu for cm. energies below the Z". So far, no 
such measurements are available. 

Our procedure to determine N^|- (J, k±^Lu) in eqs. (p!5D and (p!6| ) using the analytic expression, 
eq. ([I3|), is as follows. We begin with the compilation of the mean charged particle multiplicity 
of e+e~ hadronic annihilations, Ng+g_ versus Ec.m.i given in []T2[. These data are shown as a 
function of Ec.m. by the open symbols in Fig. ^. (For simplicity, not all the LEP-averaged results 
compiled in [|l^ are displayed in Fig. ||; however, these data are included in the procedure 
described below.) We correct these results for the contributions of c and b events so they 
correspond to our treatment of N^^^ as explained in Sect. H, i.e. to uds events only. This 
correction is a bin-by-bin multiplicative factor determined from Jetset, given by the ratio of 
Ng^g_ between uds and udscb events as a function of Ec.m.- Jetset has been found to provide 
an accurate description of the multiplicity of both the uds flavor and flavor inclusive samples in 
e"'"e~ hadronic annihilations (see for example 0, |T0|, |3^), justifying its use for this correction. 



The predictions of Jetset for Ng+g_ in udscb and uds flavor events are shown by the dashed 
and dash-dotted curves in Fig. ^. The correction factors have values between 0.92 and 0.96. 



11 



9 3J 



30 



25 - 



20 - 



15 - 



10 - 



n I I r 



10 



Unbiased qq events 



Corrected data (c & b 
quark contributions 
removed) 




^ ARGUS • 

iVjade 

□ TASSO 
A HRS 
AMY 

O LEP-averaged 



OPAL 



3NLO fit 

Jetset 7.4 udscb events 
Jetset 7.4 uds events 



Ee.„. (GeV) 



10 



Figure 3: Measurements (see [T^ for the original references) of the inclusive charged particle 
multiplicity in e"'"e~ annihilations (open symbols) as a function of -Ec.m.- The LEP results are 
averages over the LEP experiments (see |]T2[). The solid points show the corresponding results 
after corrections to remove the contributions of c and b events. For simplicity the uncertainties 
of the corrected data are not shown: these are essentially the same as those of the corresponding 
uncorrected data points. The solid curve is a fit of the next-to-next-to-next-to-leading order 
(3NL0) expression for quark jet multiplicity |T^, ^ to the corrected data. The predictions of 
the Jetset event generator for udscb and uds fiavor events are also shown. 



The corrected results for Ng^g_, corresponding to the terms Ngg(J') in eq. ( pTsD and which we 
denote Ng^'(L) in the following, are shown by the solid points in Fig. ^. We fit these corrected 
data using the next-to-next-to-next-to-leading order (3NL0) analytic expression for the scale 
dependence of quark jet multiplicity Ng [14, 36] assuming n/ = 5. Two parameters are fitted: 
an effective QCD scale parameter, Ae//., and an overall normalization constant, K. The fit 
results are given in the first row of Table |[ The result of the fit is shown by the solid curve in 
Fig. 1^. Using this fitted curve and its derivative, we numerically evaluate eq. ( |T^ ) as a function 
of 9i to determine the biased quark terms N^^- (J, k± Ly) in eqs. ([T5|) and ([T6|) . The derivative 
is evaluated by differentiating the analytic expression for with respect to y = L/2 and using 
the fitted parameter values for Ng^ (L) from Table ^ 
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Figure 4: The mean charged particle multiphcity of biased two-jet uds flavor events from Z° 
decays selected using the Durham jet finder. The results are presented as a function of p±^lu, 
related to the resolution parameter ycut of the jet finder by p_l^lu = E^s. ^/Vcut with E^is^ the 
event energy. The data have been corrected for initial state photon radiation, detector response, 
and misidentification of uds events. The uncertainties are the statistical and systematic terms 
added in quadrature. The statistical uncertainties are smaller than the symbol sizes. The solid 

The 



curve shows the prediction |T^ of the next-to-leading order (NLO) expression, eq. (|T3[) 
prediction is seem to agree with the data up to scales of about 40 GeV. The dashed curve is 
the NLO result found by altering the choice of data points used in the fit of Ng+g_ versus scale. 



6 Results 



6.1 



The biased quark jet multiplicities N^^' (J, fcj_,Ltt) 



In this section we present our direct measurement of the biased quark jet terms N^^- (L, k^^^^) 
in eq. (|T5|) . We then use these data to test the analytic expression, eq. ([T3|) . 

We apply the Durham jet finder to the sample of uds flavor events described in Sect, fusing 
a fixed value of the resolution parameter y^ut (see below). Events are retained if exactly two 
jets are reconstructed. The solid points in Fig. § show the charged particle multiplicity of the 
selected two-jet events for different choices of ycut- We relate ycut to the transverse momentum 
V^,Lu through p±^Lu = Euis. y/ycut [im, with Eyis^ the sum of the particle energy in the event. 



The data are corrected for initial state radiation, the effects of the detector, and uds event 
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Table 3: Measurements of the mean charged particle multiplicity of biased two-jet uds flavor 
events from decays as a function of the transverse momentum cutoff p±^Lu used to sepa- 
rate two- and three-jet events. The jets are defined using the Durham jet-finder. The first 
uncertainty is statistical and the second systematic. These data are shown in Fig. 



misidentification, and are plotted as a function of p±^Lu- The solid points in Fig. ^ therefore 
represent the direct measurements of the biased quark jet terms N^^- (L, k±^Lu) mentioned in 
Sect. ^. The vertical lines show the total uncertainties, with statistical and systematic terms 
added in quadrature. The evaluation of systematic uncertainties is discussed in Sect. 0. The 
statistical uncertainties are negligible compared to the systematic terms. These data are listed 
in Table |^. As part of our evaluation of systematic uncertainties, results analogous to those 
shown in Fig. ^ are determined using the Cambridge jet finder. For the Cambridge algorithm, 
the resolution parameter i/cut is related to p±^lu by the expression given above. We do not 
determine analogous results using the Lucius jet finder because it is not obvious how to relate 
the resolution parameter djoin to p±^Lu in this case. 

The solid curve in Fig. ^ shows the prediction for N^^- {L,k±^iu) found using eq. (p!3D , 
evaluated as explained in Sect. |^. The curve is seen to describe the direct measurements fairly 
well. In particular, the shape of the analytic curve is generally similar to the data, at least for 
values of scale below about 40 GeV. This suggests that the NLO result, eq. (^), is adequate 



for the purposes of our study. We return to this question in Sect |0 



The difference between the results found using the Durham and Cambridge jet finders is 
the dominant source of systematic uncertainty for the data in Fig. ^ (Sect. |^. The results 
from the Durham algorithm (standard analysis) lie closer to the NLO curve than those from 
the Cambridge algorithm, especially for small values of 6i. 

As a systematic check of the NLO prediction for N^^' {L,k±^Lu), we re-evaluated the pre- 
diction after excluding the ARGUS and 91 GeV LEP measurements from the fit of Ng+g_ 
versus -Ec.m. (see Fig. More details are given in Sect. |^. This alternate result is shown by 
the dashed curve in Fig. The difference between the solid and dashed curves is seen to be 
similar to the difference between the solid curve and symbols. This suggests that the agree- 
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ment of the analytic expression eq. ([T3|) with the direct measurements is quite reasonable, once 
sources of systematic uncertainty are considered. 



6.2 Unbiased gluon jet multiplicity versus scale 



Fig. H shows our results for the unbiased charged particle multiplicities of gg events, N^^-. The 
solid points in Fig. ^ are obtained from eq. (|TB|) using the direct measurements of N^^- (L, k^^^u), 



i.e. the solid points in Fig. The small horizontal bars indicate the statistical uncertainties. 
The vertical lines show the total uncertainties, with statistical and systematic terms added in 
quadrature. The evaluation of systematic uncertainties is discussed in Sect. |^. The asterisks 
and open symbols in Fig. Pd show the corresponding results from eqs. ( |l5l) and (|T^) using the 
calculated expressions for N^|- (J, k±^Lu) from eq. (|I3]). Thus the asterisks in Fig. ^ are derived 
using the solid curve in Fig. H, for example. The two sets of results from eq. ([T5|) (solid points 
in Fig. I^a and asterisks in Fig. ^) are seen to be very similar. This is consistent with the 



observation in Sect. 6.1 that eq. ([T3|) is generally adequate for the purposes of our study. 



The results from eq. (|T^) are plotted as a function of Q = p±^lu and those from eq. (|TBp as 
a function of Q=p±,Le- The variable p±^Lu extends over a smaller range of scale than p±^Le 
because the factor of s in the denominator of the Lund definition of transverse momentum 
is always larger in Y events than the factor of Sqg in the Leningrad definition, see eqs. (|Tl|) 
and (p!2|). This effect is also visible in Fig. §3 and Table |l|, i.e. p±^Lu has a maximum of 
30.5 GeV compared to 53.1 GeV for p±^Le- Note that with 6i > 35° as in our standard analysis, 
the lower bound on the gluon jet scale is 11.1 GeV for p_l^lu and 16.8 GeV for p_L^Le, see Fig. ||b. 

Included in Fig. ^ are direct measurements of the inclusive charged particle multiplicity of 
unbiased gluon jets from the CLEO and OPAL Collaborations, currently the only results which 
utilize a hemisphere definition of gluon jets. These results are shown by the triangular symbols. 
The three data points at scales of about 4-7 GeV are derived from the hadronic component of 
T(1S) — > 'ygg events ^j. The scale Q for these data is given by the invariant mass of the hadronic 
system. Similarly, T(3S) — > 7x^2(2?) — > 'jgg events provide the result at Q = 10.3 GeV 0, with 
the scale given by the Xb2{'^P) niass. The measurement at Q = 80.2 GeV |jll| is determined using 



hadronic Z" decays: Z° qqgind., where gind. refers to a gluon jet hemisphere recoiling against 



two almost coUinear identified quark jets q and q in the opposite hemisphere |33|. For our 



purposes the g^d. measurement from [|r^ has been multiplied by a factor of two both for the 
multiphcity and energy scale so it corresponds to gg "two-jet events" analogous to the other 
data in Fig. ^. 



Fig. ^ also shows the predictions of the Herwig p8| , ^ and Jetset Monte Carlo event 



generators for the inclusive charged particle multiplicity of gg events versus Q = Ec.m.- For these 
results we use version 6.2 of Herwig with the following changes to the default parameter 



set: the parameter specifying the angular smearing of udsc quark clusters, CLSMR(l), is set 
to 0.40, the parameter controlling the mass spectrum of b quark clusters, PSPLT(2), is set 
to 0.33, and the parameter controlling the production of decuplet baryons, DECWT, is set 



^We parametrize the data in Fig. ^ using a polynomial. We then evaluate this polynomial at the scale values 
corresponding to our bins in 9i (see Table |^) to determine the gluon jet multiplicities using eq. ([l5|). 
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Figure 5: 



The mean charged particle multiphcity of unbiased gg events as a function of scale, 
(a) Results from eq. (p!5| ) using the measured biased quark jet terms N^^- (L, The small 

horizontal bars indicate the statistical uncertainties. The vertical lines show the total uncer- 
tainties, with statistical and systematic terms added in quadrature. The sohd curve is a fit 
of the 3NL0 expression for gluon jet multiplicity [|14[ ^ to the data, (b) The correspond- 
ing results from eqs. (|T^ and (|TB|) using the NLO expression for the biased quark jet terms 
Ng^- (J, k±^Lu), eq. (0), where J = L for the N^^' {k±^Lu) results and J = Lqq for the N^^' {k±^Le) 
results. For both (a) and (b), the triangles show measurements of the inclusive charged particle 
multiplicity of unbiased gluon jets from the CLEO ^ and OPAL Collaborations. The 
predictions of the Herwig and Jetset Monte Carlo event generators are also shown. 



16 



to 0.70. Herwig 6.2 with these changes is found to provide a good description of the global 
features of the OPAL data including those of flavor selected samples. 



From Fig. ^ it is seen that the Monte Carlo predictions describe the direct measurements of 
unbiased gluon jet multiplicity (triangle symbols) well. The results from the present analysis 
based on eq. (p!5| ) (solid points in Fig. ||a and asterisks in Fig. ||b) are also well described by 
the Monte Carlo curves. In contrast, the results from eq. (0) (open symbols in Fig. ^d) are 
generally well above the Monte Carlo predictions, and - if extrapolated to lower and higher 
energies - appear inconsistent with the direct measurements from CLEO and OPAL as well. 
We note that both Jetset and Herwig have been found to provide a good description of gluon 
jet properties, including multiplicity, in other studies (see for example 0)- We 



therefore conclude that the equation based on the Lund definition of the gluon jet scale, eq. (^, 
yields results which are more consistent with other studies than the equation based on the 
Leningrad definition, eq. (^). Our data thus allow a fairly clear discrimination between the two 
possibilities eqs. (||) and (^. We henceforth restrict our analysis of gluon jets to the former 
set of results. Numerical values for the gluon jet measurements based on eq. (H) (and (|15D) 
are presented in Table ^ The analogous results based on eq. (^) (and (0)) are presented in 
Table For information, we include results corresponding to all 6i bins in these tables, not 
just the bins above 35° used in our standard analysis. Note that the theoretical formalism is 
not expected to be valid for small values of 6i. This explains the negative values of multiplicity 
determined for the first bin, < < 5°. 



6.3 The ratios r, r^^^ and r^^^ 
6.3.1 Determination of the ratios 

Having established that the formalism of [T^ yields consistent results with the direct measure- 



ments of unbiased gluon jet multiplicity if the equation for three-jet event multiplicity based 
on the Lund definition of transverse momentum is employed (Sect. |6.2|) , we proceed to a com- 
parison of the unbiased gluon and quark jet terms and to a test of the theoretical expressions 
for the ratios r, r^^^ and r^'^\ see eqs. (P-(^. For this purpose we use the gluon jet data 
shown by the solid points in Fig. |^a since these results are derived using the direct measure- 
ments of N^l' (L, k±^Lu) rather than the NLO expression, eq. (pIS]) . They therefore rely on fewer 
theoretical assumptions. 

We employ the following procedure to determine r, r^^^ and r*^^-* from experiment. The 
corrected unbiased quark and gluon jet multiplicities, shown by the solid points in Figs. |] 
and are separately fitted using the 3NL0 expressions for N, and N^, [|I^, respectively. 
To better constrain the results, the measurements of inclusive gluon jet multiplicity from CLEO 
at 10.3 GeV and OPAL at 80.2 GeV are included in the gluon jet fit. The fitted parameters are 
Ag//. and the normalization constant K. (For the quark jet data, this is the same fit discussed 
in Sect. |^.) The results of the fits are shown by the solid curves in Figs. |^ and ^ja. The fitted 
parameters and corresponding per degree-of-freedom are listed in Table |^. The per d.o.f. 
and overall description of the data are seen to be reasonable in both cases. The ratio of the 
fitted expressions for N^, and Ng defines r, eq. ([|). We calculate the first and second derivatives 
of the analytic equations for Ng and N^ with respect to y and evaluate the resulting expressions 
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Q (GeV) 


6*1 (degrees) 




r 


1.23 


0-5 


-8.8 ±0.3 ±6.4 


-3.1 ±0.1 ±2.1 


2.10 


5-10 


0.9 ±0.1 ±4.8 


0.3 ±0.0 ±1.5 


3.29 


10-15 


4.6 ±0.1 ±2.4 


1.00 ±0.01 ±0.56 


4.88 


15-20 


6.9 ±0.1 ±1.8 


1.25 ±0.01 ±0.34 


6.47 


20-25 


8.0 ±0.1 ± 1.9 


1.28 ±0.02 ±0.31 


8.01 


25-30 


8.9 ± 0.2 ± 1.9 


1.28 ±0.02 ±0.27 


9.52 


30-35 


9.7 ± 0.2 ± 1.9 


1.29 ±0.02 ±0.26 


11.1 


35-40 


10.6 ±0.2 ±1.8 


1.31 ±0.03 ±0.22 


12.6 


40-45 


11.0 ±0.3 ±1.7 


1.29 ±0.03 ±0.20 


14.1 


45-50 


11.2 ±0.3 ±1.7 


1.25 ±0.03 ±0.19 


15.5 


50-55 


13.0 ±0.3 ±1.9 


1.39 ±0.04 ±0.20 


17.0 


55-60 


13.2 ±0.4 ±1.8 


1.35 ±0.04 ±0.19 


18.4 


60-65 


13.3 ±0.4 ±1.7 


1.31 ±0.04 ±0.17 


19.9 


65-70 


14.4 ±0.5 ±1.2 


1.38 ±0.05 ±0.12 


91 A 




1 A c, _|_ n c; -1- 1 1 




22.9 


75-80 


15.00 ±0.52 ±0.87 


1.353 ±0.047 ±0.079 


24.0 


80-85 


15.53 ±0.57 ±0.94 


1.371 ±0.050 ±0.083 


25.4 


85-90 


15.8 ±0.6 ±1.5 


1.36 ±0.05 ±0.13 


26.7 


90-95 


15.9 ±0.6 ±1.9 


1.34 ±0.05 ±0.16 


27.9 


95-100 


16.5 ±0.7 ±2.5 


1.37± 0.06 ±0.21 


28.7 


100-105 


17.9 ±0.7 ±2.0 


1.46 ±0.06 ±0.16 


29.6 


105-110 


17.6 ±0.7 ±2.3 


1.42 ±0.06 ±0.18 


30.2 


110-115 


19.1 ±0.7 ±1.6 


1.53 ±0.06 ±0.13 


30.5 


115-120 


18.2 ± 0.7 ± 1.9 


1.45 ±0.06 ±0.15 



Table 4: Measurements of unbiased gluon jet multiplicity N^^ as a function of energy scale 
Q =p±^Lu, obtained using eq. ([T5|) . The corresponding bins of 9i in Y events are also indicated. 
The results for 6i > 35°, listed in the bottom part of the table, are used in our standard analysis. 
These data are shown by the solid points in Fig. For information, we also list the results 
obtained for smaller values of 6i. Note that the theoretical formalism used to determine these 
results is not expected to be valid for small 6i. This explains the negative value of multiplicity 
derived for 9i < 5°. The corresponding results for the multiplicity ratio r between gluon and 
quark jets are also given. The results for r for 6i > 35° are shown by the solid points in Figs. ^ 
and l^a. 
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Q (GeV) 


6i (degrees) 




1.59 


0-5 


-1.2 ±0.3 ±6.5 


2.90 


5-10 


6.8 ±0.1 ±4.8 


4.66 


10-15 


9.5 ±0.1 ±2.1 


7.10 


15-20 


11.6 ±0.1 ± 1.5 


9.58 


20-25 


12.9 ±0.1 ± 1.5 


12.0 


25-30 


14.1 ± 0.2 ± 1.3 


14.3 


30-35 


15.2 ± 0.2 ± 1.2 


16.8 


35-40 


16.4 ± 0.2 ± 1.2 


19.2 


40-45 


17.1 ±0.3 ± 1.3 


21.8 


45-50 


17.6 ± 0.3 ± 1.4 


24.0 


50-55 


19.6 ± 0.3 ± 1.7 


26.5 


55-60 


20.0 ± 0.4 ± 1.5 


29.0 


60-65 


20 2±04±14 


31.7 


65-70 


21.44 ±0.47 ±0.99 


34.6 


70-75 


21.65 ±0.49 ±0.99 


37.4 


75-80 


22.25 ±0.52 ±0.95 


39.2 


80-85 


22.86 ±0.57 ±0.96 


42.0 


85-90 


23.3 ± 0.6 ± 1.5 


44.8 


90-95 


23.4 ± 0.6 ± 1.8 


47.6 


95-100 


24.1 ±0.7 ±2.5 


49.0 


100-105 


25.6 ± 0.7 ± 1.9 


51.1 


105-110 


25.4 ±0.7 ±2.1 


52.6 


110-115 


26.9 ± 0.7 ± 1.6 


53.1 


115-120 


26.0 ± 0.7 ± 1.9 



Table 5: Measurements of unbiased gluon jet multiplicity N^^' as a function of energy scale 
Q =p±^Le, obtained using eq. (|T6|). Note that the theoretical formalism used to determine these 
results is not expected to be valid for small 6i. This explains the negative value of multiplicity 
derived for < 5°. The results corresponding to 6i > 35° are shown by the open points in 
Fig. i). 
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Figure 6: Results for the ratio r of the mean charged particle multiplicities between unbiased 
gluon and uds quark jets as a function of scale. The central result is shown by the solid curve; the 
shaded and open bands show the corresponding statistical and total uncertainties, respectively 
(the statistical uncertainties are barely visible). The short-dashed curve is an extrapolation of 
the fit. The solid points show measurements based on gluon jet results from the present analysis. 
The triangle and square symbols show results based on direct measurements of multiplicity in 
unbiased jets from the CLEO % §, ARGUS Q and OPAL Collaborations. The dash- 



dotted and long-dashed curves show theoretical predictions based on analytic ||T^, |T5[ and 
numerical [O techniques. 



using the corresponding parameter values from Table 0. The ratios of these terms define r^^^ 
and r*^^), eqs. (0) and @. 

Our results for r are shown in Fig. and those for r^^^ and r*^^^ in Fig. |^. The central 
results are indicated by solid curves. The shaded bands show the statistical uncertainties. The 
overall uncertainties, with statistical and systematic terms added in quadrature, are shown by 
the open bands. The bands cover the energy range of data in the gluon jet fit, from 10.3 to 
80.2 GeV. The short-dashed curves show extrapolations of the results to lower and higher scales. 
At 30 GeV, a typical scale for the gluon jet data in these figures, the ratios are determined to 
be: 

r (30 GeV) = 1.422 ± 0.006 ± 0.051 , (17) 
r(^)(30GeV) = 1.761 ± 0.013 ± 0.070 , (18) 
r(2)(30GeV) = 1.98 ± 0.02 ± 0.13 , (19) 

where the first uncertainty is statistical and the second systematic. At 80 GeV, corresponding 
to the highest energy gluon jet data in the fits, the results are: 

r (80 GeV) = 1.548 ± 0.008 ± 0.041 , (20) 
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Figure 7: (a) and (b): The ratios of slopes r*^^-* and of curvatures r*^^^ between unbiased gluon 
and uds quark jets as a function of scale. The curves and bands have the same meaning as in 
Fig.§ 
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rW(80GeV) = 1.834 ± 0.016 ± 0.088 , (21) 
r(2)(80GeV) = 2.04 ± 0.02 ± 0.14 . (22) 

Analogous results for r versus scale were previously published in [18], and for r^^^ in [jl^, [T9[|. 



These earlier studies were based on biased jet samples, however. Therefore the quantitative 
results from those studies cannot be compared to theoretical predictions without ambiguity. 
Our result for r^"^^ is the first experimental determination of that quantity. 

Included in Fig. ^ are the results for r derived from the individual gluon jet measurements 
in Fig. The solid points in Fig. ^ are obtained by dividing the results shown by the solid 
points in Fig. ^ by the results of the 3NL0 parametrization of quark jet multiplicity at the 
corresponding scales. Numerical values of these results for r are compiled in Table The 
triangle and square symbols in Fig. ^ are obtained as follows: (i) for the OPAL data point 



at 80.2 GeV, the published result ||rT| is used, based on uds quark events analogous to the 
present study; (ii) for the three CLEO points at 4-7 GeV, the published results [|] are used 
after applying a correction obtained from Jetset to remove the contributions of c events from 
the quark jet samples; (iii) for the point labelled "CLEO/ARGUS," the CLEO measurement 
of gluon jet multiplicity at 10.3 GeV ^ is divided by the ARGUS measurement of N^^g. at 
10.5 GeV 1^ (see Fig. ||), after applying a correction to account for c and b events in the quark 



jet sample and for the small energy difference between the gluon and quark jet results. 

As a consistency check, we also determined the ratios r and r^^^ using the technique employed 



m 



T9| . The corrected quark jet results (solid points in Fig. ^ were fitted using a polynomial 
expression: 

Ngg = ao + ai L + a2 L'^ , (23) 

with L given in eq. @, and where ao, ai and 02 are fitted parameters. We find ao = 13.20±0.67, 
ai = — 2.91±0.12 and a2 = 0.2818±0.0057, where the uncertainties are statistical. The (d.o.f.) 
of the fit is 3.74 (15). We then fitted the 19 gluon jet measurements in our analysis (solid points 
in Fig. ^ and the 10.3 and 80.2 GeV results from CLEO and OPAL) using the expression: 

= bo + h , (24) 

with the results 60 = -4.07 ± 0.27, h = 1.759 ± 0.030 and (d.o.f.) = 31.5 (17). The result 



for r, given by the ratio of eqs. (p^ and (23), is indistinguishable from the solid curve in Fig. |g. 
The result for r^^\ obtained by differentiation of eq. (|2^), is, by construction, a constant: 
r^^^ = bi = 1.759 ± 0.030 (stat.). This result is consistent with that shown in Fig. |^, see also 
eq. ([T8|) . Not only is r*^^^ constant with this method, but r^^^ =bi =r^^\ i.e. no distinction is 
made between the ratios of slopes and curvatures. Therefore this alternate technique does not 
provide as much information as our standard procedure described above. 

To obtain a consistency check of our result for r^^\ we applied the following procedure. 
The 19 gluon jet measurements were clustered into 16 groups of four contiguous measurements. 
Thus, labelling the data points from 1 to 19 in order of increasing Q, the first group is composed 
of data points 1-4, the second of 2-5, etc. A straight line fit was performed for each group: 

(N,g), = (a,),L+(6,), z = l,16 , (25) 

with ttg and bg the fitted parameters. By fitting to four points at a time, rather than e.g. to two 
points only, the influence of bin-to-bin fluctuations is reduced. We then evaluated the fitted 
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polynomial expression eq. (p3|) at the scale values of the gluon jet measurements to obtain 19 
corresponding results for quark jets, which were then fitted in the same manner: 

{^gg)i = ia,)iL + ibg), 1 = 1,16 . (26) 

A scale value was associated with each slope point {ag)i or (ag)j using the arithmetic mean 
of the four Q values in the fits. As a systematic check, we also defined these scales using the 
weighted means of the Q values of the four associated gluon jet measurements. This alternate 
method yielded consistent results with those presented below. The slopes {ag)i and (aq)j were 
in turn fitted using straight lines: 

ag=CgL + dg (27) 
ttg = Cg L + dg (28) 

with Q and di [i = q, g) the fitted terms. The ratio r^^-* found using this technique is a constant 
given by r^"^^ = Cg/cg. We find r*^^) = 2.25 ± 0.16 (stat.), consistent with our result in Fig. |^, see 
also eq. ([T9|) . 



The procedure described in the previous paragraph can also be used to obtain a second 
consistency check for r^^\ Taking the weighted mean of the 16 results for the ratio of slopes 
rl^^ = {ag/aq)i yields r^^^ = 1.68 ± 0.02 (stat.), consistent with the result in eq. ([TsD to within 
about 1.0 standard deviation of the total uncertainties. 



We note that the method to determine r^^^ and r*^^^ based on eqs. (p5|)-(p8|) is more direct 
than our standard procedure because no functional form is assumed for how the gluon jet 
multiplicity varies with scale. The resulting uncertainties are considerably larger, however, and 
no information on the scale dependence of r*^^^ and r*^^^ is obtained. Therefore we retain the 
analysis based on the 3NL0 parametrizations of multiplicities for our standard results. 

Analytic expressions for r, r 

(1) and r(2), valid to 3NL0, were recently presented by Capella 
et al. fl^. All three ratios are predicted to converge to the ratio of QCD color factors, Ca/Cf = 
2.25, in the asymptotic limit of large jet energies. This convergence is predicted to be more 
rapid for r^^^ than r |14, 18| and yet more rapid for r*^^^ Our data are in agreement with 



this prediction, i.e. r < r*^^^ < r^"^^ < 2.25 for the scales accessible in our study, see eqs. ([T7|)-(P^. 
Note that we also observe the hierarchy r < r'-^-' < r*^^^ using the alternate analysis strategies 
summarized by eqs. (H)-(H). 

To make the comparison of our results for r, r^^-* and r*^^^ more clear, the data from Figs. |] 
and 1^ are displayed together in Fig. The dominant uncertainties in our analysis are system- 
atic (see Sect. ^ and are correlated between the three ratios. Therefore, in Fig. Pd, we present 
results for the differences r*^^^ — r, r^^^ — r and r^^^ — r^^\ for which correlated systematic uncer- 
tainties partially cancel. At 30 GeV, we find: 



^(2) _ ^(1) 



(30 GeV) = 0.56 ±0.02 ±0.14 , (29) 
(30 GeV) = 0.341 ±0.014 ±0.075 , (30) 
(30 GeV) = 0.221 ±0.024 ±0.070 , (31) 



where the first uncertainty is statistical and the second systematic. These results demonstrate 
that r*^^) exceeds r^^^ by 3.0 standard deviations of the total uncertainties, with yet more 
significant deviations from zero observed for the other two differences of ratios. 
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Figure 8: (a) The experimental results for r, r^^^ and r^"^^ from Figs, ^and 0, collected together; 
(b) the differences between ratios r^"^^ — r, r^^^ — r and r*^^) — r^^^ as a function of scale. The 
central results are shown by solid curves; the shaded and open bands show the corresponding 
statistical and total uncertainties, respectively. 
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6.3.2 Comparison to theoretical expressions 



The 3NL0 analytic predictions of Capella et al. [0] for r, r*^^^ and r^^^ are shown by the long- 
dashed curves in Figs. || and |^. For Q = 30 GeV, the ratios are predicted to be 1.73, 1.87 
and 2.00, respectively, assuming nf = 5 and A = 0.20 GeV. Thus, the analytic predictions for r 
and r^^) exceed the corresponding experimental results in eqs. (|r^)-([T8|) by about 22% and 6%, 
while the theory agrees with the data, eq. (|1^), for r^'^\ Therefore the analytic predictions 
agree better with the data the higher the level of differentiation of the multiplicity curves. This 
suggests that higher order corrections are smaller for r^^^ than r^^\ and for r^^^ than r. 

We note there are ambiguities in the analytic calculation related to the effective number of 
active quark flavors and the scale Q at which to evaluate the strong coupling strength. The 
3NL0 predictions shown here are based on n/ = 5 (as stated above) and Q = Ec.m.- Other 
plausible choices such as ^^/ = 4 and Q = Ec.m./^ change the analytic prediction for r by as 
much as about 10% and can bring the theoretical curve into better agreement with the data, 
see Fig. 2 of for example. 



A second QCD prediction for r versus scale was recently presented by Lupia and Ochs |13]. 



This prediction is based on numerical, rather than analytic, techniques, incorporating more 
accurate phase space limits for the emission of soft gluons and a more complete treatment of 
energy conservation than the analytic expression. Unlike the analytic result, a prediction is 
presented only for r, however, not for r*^^) or r^'^\ We derive a prediction for r^^^ from the 
numerical calculation in the following manner. We begin with numerical predictions for the 
individual gluon and quark jet multiplicities,^ given in steps of 0.02 units in y with y = In (Q/Qc) 
where Qc = 0.507 GeV [|13[. We make linear interpolations between the points and use the 
slopes of these lines to estimate dNg/dy and dNg/dy, which are then used to determine r^^^ as 
a function of scale Q. We attempted to extend this technique to obtain a prediction for r^'^^ 
as well. The derived results for r^"^^ were very scattered from point to point as a function of Q 
because of limited numerical precision. We therefore do not quote a prediction for r^"^^ based 



on the calculation of [|T3| 



The predictions of the numerical calculation for r and r^^^ are shown by the short-dashed- 
dotted curves in Figs. ^ and 0a. At 30 GeV, the numerical calculation predicts r to be 1.45, 
in agreement with our measurement at that scale, eq. (p^Tf ). The corresponding result for r*^^^ 
is 1.64, about 7% below the data, eq. (|18|). For r, the agreement with data is therefore better 
for the numerical method than for the analytic calculation. The level of agreement of the two 
predictions with data is similar for r^^\ The difference between the analytic and numerical 
results for r (long-dashed and short-dash-dotted curves in Fig. |^) is believed to arise from the 
differences in phase space limits and energy conservation noted above. For further discussion, 
see 
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Finally, theoretical predictions for r, r^^^ and r*^^^ can be derived from the formalism of Eden 
and Gustafson [|T5|, using eq. (p!^). To obtain a prediction for r^^\ we begin with the experi- 



mental results for the slope of quark jets, dNqg(L)/dL, determined in the manner described in 
Sect. |6.3.1| for dNqq{y)/dy, with y = L/2. We insert this result in eq. (|1^ to obtain a predic- 



tion for d'Ngg{L)/dL. Dividing this result by the measured d'Ngg{L)/dL yields the prediction 
^We thank Wolfgang Ochs for providing these resuhs. 
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for r^^\ To obtain a prediction for r^'^\ we differentiate eq. (T|) with respect to L and follow an 
analogous procedure to that described for r^^\ To obtain a prediction for r, we fix the gluon jet 



multiplicity at 80.2 GeV to the gind. measurement from OPAL |TT[. We then use the prediction 
for dNgg{L)/dL determined as described above to evaluate the gluon jet multiplicity at other 
scales and divide the result by the 3NL0 expression for the quark jet multiplicity fitted to data 
(solid curve in Fig. |^). 

The predictions obtained in this manner are shown by the long-dash-dotted curves in Figs. |] 
and 1^. The results are seen to be in good overall agreement with the data. We note, however, 
that these predictions are based on the experimental measurements of quark jet multiplicities. 
Furthermore, the prediction for r utilizes the measured gluon jet multiplicity at Q = 80.2 GeV 
as stated above. Therefore, the predictions we derive based on are not independent of the 
data. 



6.4 Measurement of Ca/Cf 

Measurements of the energy dependence of unbiased particle multiplicities in gluon and quark 
jets can also be used to determine an effective value of the ratio of QCD color factors Ca/Cf. 



Such a possibility was first discussed in |]T8|. To perform this measurement, we integrate eq. (|T^) 
to obtain: 

N,,(L)=N,,(Lo) + ^^(L) , (32) 
with Ngg (Lq) an integration constant which serves as the reference point for N^^ (L), and where 

aoCr\ dNqg{x) 



fL — Cg+Cq 

J Lo — Cg+Cg 



X J dx 



dx . (33) 



We evaluate eq. ( |5B| ) numerically using the experimental results for dN^g (Q) / dQ discussed 
in Sect. |6.3.1| . As the reference point N^^ (Lq) in eq. (|32D, we use the direct measurement 
of gluon jet multiplicity from OPAL 0. Thus Lq = 2 In (Qo/A), with go = 80.2 GeV and 
'Ngg (Lq) = 28.56 |11|]. For A we used the fitted result for Ag//. for quark jets from Table |^, 
A = 0.19 GeV. Note that the results for Ca/Cf are independent of this choice. 

A one parameter fit of eq. ( |52D is made to the gluon jet measurements in Fig. |^a, excluding 
the CLEO data below 7 GeV, to obtain 

C 

= 2.232 ±0.008 (stat.) ±0.108 (syst.) ±0.082 (ref.) , (34) 

Cp 

where the first uncertainty is statistical. The second uncertainty is a systematic term evaluated 
as described in Sect. The third uncertainty is an additional systematic term defined by the 
difference observed if the CLEO measurement of Ngg at 10.3 GeV is used as the reference point, 
rather than the OPAL measurement at 80.2 GeV. The (d.o.f.) of the fit is 29.5 (17). The 
result of our fit, eq. (p5), is consistent with the QCD value Ca/Cf = 2.25. 



The fitted expression of eq. ( p2D yields a curve which is almost identical to the 3NL0 curve 
in Fig. I^a: therefore we do not display the fitted curve of eq. (^) in addition. In Fig. ^ 
we present the result for the event multiplicity in Y events. The data in this figure are the 
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Figure 9: The mean charged particle multiphcity of three-jet uds flavor Y events from Z° 
decays compared to the analytic expression, eq. (^), for which the gluon jet multiplicity term 
{k±^Lu) is determined from a one parameter fit of eq. (0) to unbiased gluon jet data and 
the quark jet term N^g (L, k±^Lu) is measured. The open band shows the evaluated uncertainty 
of the analytic result. Note the theory is not expected to be applicable for small values of 9i 
because the events in this case do not have a well defined three jet event structure. 



same as in Fig. ^a. The curve is obtained from eq. (^), with N^^ (fc±,Lu) the fitted expression 
of eq. ( p2D and N^^ (L, the measured result from Sect. The open band shows the 

uncertainty of the analytic result, evaluated using the uncertainty from eq. ( |5^ ) in addition to 
those shown for the Ng^ (L, terms in Fig. ^. The curve is seen to provide an accurate 

description of the data to within its uncertainty for 6i values larger than approximately 10°. 
This represents a considerable improvement compared to the results of based on earlier 



theoretical formalism [|T7l, in which the corresponding fitted analytic expression was found to 
describe the data for 6i values larger than about 60° only (see Fig. 4 of 



To test our technique of determining Ca/C^, we apply our analysis to parton level Monte 
Carlo events generated with a cm. energy of 10 TeV. The parton level is defined using final- 
state partons, i.e. those which are present after the termination of the parton shower. We 
use Herwig version 6.2 with the parameters discussed in Sect. p.'J\ . Applying our analysis to 
parton level events with a very large energy tests our method to determine Ca/Cf since the 
parton level asymptotic result must be consistent with 2.25 for the analysis to be considered 
valid, see for further discussion. The parton level Herwig events at 10 TeV yield Ca/Cf = 



2.285 ± 0.010 (stat.) ± 0.047 (ref.), where the central value is defined using the Herwig parton 
level prediction for N^^ at 80 GeV as the reference point. The first uncertainty is statistical while 
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the second uncertainty is defined by the difference found if the Herwig prediction for N^^ at 
10 GeV is used as the reference point instead (i.e. this is analogous to how the third uncertainty 
of eq. ( p^ ) is defined). Thus our one parameter fit method yields an asymptotic result of 
2.25 to within the uncertainties associated with the analysis procedure, implying that the 
similarity of our measurement eq. (0) with the QCD value of Ca/Cf is not merely coincidental. 
For information, the parton level Herwig result at -Ec.m. =91 GeV is 2.03 ± 0.01 (stat.) if the 
prediction for N^^ at 80 GeV is used as the reference point and 2.35±0.02 (stat.) if the prediction 
at 10 GeV is used. The analogous results for hadron level Herwig events are 2.34 ± 0.01 and 
2.07 ±0.02, respectively. We therefore observe a much stronger dependence on the choice of the 
normalization point for the 91 GeV Monte Carlo results than for the 10 TeV results or for the 
data. The Herwig results at 91 GeV are nonetheless consistent with 2.25, once the systematic 
uncertainty related to choice of the normalization point is considered. 



7 Systematic uncertainties 

To evaluate systematic uncertainties, the analysis was repeated with the following changes to 
the standard analysis. The full differences between the standard results and those found using 
each of these changes were used to define symmetric systematic uncertainties. For items 1, 5 
and 7, listed below, the largest of the described differences with respect to the standard result 
was assigned as the systematic uncertainty. The uncertainties were added in quadrature to 
define the total systematic uncertainties. 

For the measurement of the gluon jet multiplicity N^^' (Fig. the following changes were 



made. Items 1-3 and 5 in this list were also applied to the study of two-jet events, see Sect. |0. 



The items are listed in roughly decreasing order of the size of their contribution to the total 



uncertainties of the N^^' measurements. 



The jet finding was performed using the Cambridge and Lucius jet finders, rather than the 
Durham jet finder. For the Lucius jet finder, we did not perform a direct measurement 
of the biased quark jet terms ^a^' {J, k±^Lu) since it is not obvious how to relate the 



resolution scale djoin to p±^lu as mentioned in Sect. pAj Instead we used the analytic 
expression, eq. (0), for the determination of the gluon jet multiplicities N^^' (L). The 
difference between the Durham and Cambridge results was used to define the systematic 
uncertainty due to jet finding for the two jet event measurements shown in Fig. ^. 

2. Tracks selected for the uds tagging procedure were required to have a signed impact 
parameter which satisfied dsign/ad^^g^> 1-5, rather than dsign/crd^.^„>S.O; at the same 
time we required | cos(6'thrust)| to be less than 0.8, rather than 0.9. The purpose of these 
two changes was to increase the purity of tagged uds events; the estimated uds event 
purity increased to 85.1% from 78.5% while the number of selected events decreased to 
8610 from 22 365. The more restrictive condition on | cos(^thrust)| yields a higher uds 
purity because a larger fraction of tracks in an event traverses the silicon microvertex 
detector. 

3. The Herwig Monte Carlo event generator was used to determine the corrections for initial- 
state photon radiation, detector response, and uds event and gluon jet misidentification. 
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rather than Jetset. For this purpose, we used version 5.9 of the program ||3^ with the 
parameter values in [^], except that the cluster mass cutoff CLMAX was increased from 
3.40 to 3.75 GeV/c^ to improve the model's description of the mean charged particle 
multiplicity in inclusive hadronic decays. 

4. For our measurements of the gluon jet multiplicities N^^' based on eq. (]T5|), we used the 
analytic expression, eq. ([131) , rather than the direct measurement of N^^- (L, k±^Lu)- 

5. Charged tracks alone were used for the data and Monte Carlo samples with detector 
simulation, rather than charged tracks plus electromagnetic clusters. As an additional 
check on the track selection, the standard analysis (tracks and clusters) was varied by 
increasing the minimum transverse momentum of charged tracks with respect to the beam 
axis from from 0.05 GeV/c to 0.15 GeV/c. 

For the two-jet events (Fig. |^), the dominant systematic term was from item 1 at low p±^Lu and 
from item 3 at intermediate and high p±^Lu- 

In addition to the above, the following changes were made to evaluate systematic uncer- 
tainties for the ratios r, r'-^-* and r*^^-*. 

6. The ARGUS and 91 GeV LEP measurements of N^+g, were excluded from the fit of quark 
jet multiplicity (solid curve in Fig. These two data points have the smallest uncer- 
tainties of the results in the fit and were excluded to potentially maximize a systematic 
variation. The results of the alternate fit are Ae//. = 0.261 ±0.061, K = 0.154 ±0.014 and 

(d.o.f.) = 1.5 (14). Besides the evaluation of a systematic uncertainty for r, r^^^ and 
this alternate result for N^^' (L) was used to derive the NLO prediction shown by the 
dashed curve in Fig. ^j. 

7. Y events with 50° < < 120° were used in the fit of gluon jet multiplicity versus scale, 
rather than 35° < < 120°; as an alternate check, the ratios were determined after in- 
creasing the inclusive gluon jet measurement from CLEO at 10.3 GeV by its one standard 
deviation total uncertainty, then decreasing it by one standard deviation, and repeating 
the same operation for the inclusive gluon jet measurement from OPAL at 80.2 GeV. 

For the ratios, the largest systematic terms were from items 1, 6 and 7. The contributions 
of items 2 and 5 were by comparison negligible, with the other items intermediate. For our 
measurement of Ca/Cf (Sect. |6.4|) , increasing and decreasing the inclusive OPAL measurement 
at 80.2 GeV by one standard deviation (item 7) provided the largest contribution to the overall 
uncertainty. Item 1 provided the second largest contribution while the contributions of the 
other items were much smaller. 

For the data listed in Tables | and |] and the corresponding results in Figs. |, |^ and ^ja, 
the systematic uncertainty evaluated for each bin was averaged with the results from its two 
neighbors to reduce the effect of bin-to-bin fluctuations. The single neighbor was used for bins 
on the endpoints of the distributions. 
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8 Summary 



We have presented measurements of the mean charged particle multiphcity of three-jet "Y 
events" as a function of the opening angle between the two lowest energy jets, 9i, and of two- 
jet events as a function of the resolution scale separating the two- and three-jet event classes. 
The measurements were performed using the Durham jet finder. The Cambridge and Lucius jet 
finders were used as systematic checks for the three-jet measurements and the Cambridge jet 
finder for the two-jet measurements. The two- and three-jet events were selected from a sample 
of Z° decays to light flavor (u, d or s) quark-antiquark pairs produced in e^e~ annihilations. 
The restriction to light quark flavors improves the correspondence of the data to theoretical 
predictions. 



We use these data to test recent theoretical formalism |T^, |T^ for particle multiplicity in jets 



which accounts for biases introduced by the jet finding criteria. We find that the theoretical 
prediction for particle multiplicity in biased two-jet events, eq. (|13|), agrees fairly well with 
our two-jet event data. In conjunction with the results from Y events, the two-jet data are 
used to extract the multiplicity of unbiased gluon jets at a variety of scales. Of two possible 



forms proposed in ||T6| to obtain this information - one [|T5[ based on the definition of gluon 
jet transverse momentum from the Lund group and the other from the Leningrad group - 
we find clear preference for the former because of the consistency of the derived results with 
direct measurements of unbiased gluon jet multiplicity from the CLEO [0, |[ and OPAL |P|- [pl| 
Collaborations as well as with Monte Carlo expectations. 

The unbiased gluon jet multiplicities extracted with this technique, which range in scale 
from 11.1 to 30.5 GeV, are compared to corresponding results from light unbiased quark 
jets. In conjunction with direct measurements of unbiased gluon jet multiplicity at 10.3 and 
80.2 GeV from CLEO and OPAL, we determine the ratio of charged particle multiplicities 
between gluon and quark jets, r = N^/Ng, and the corresponding ratios of slopes and of cur- 
vatures, r^^) = {dNg/dy)/{dNq/dy) and r^^^ = (d^Ng/d?/^)/(d^Ng/d?/^), as a function of energy 
scale y = \n{Q/A), with Q the jet energy and A the QCD scale parameter, set to 0.20 GeV 
in our study. At 30 GeV, a typical scale in our analysis, we find r = 1.422 ± 0.006 ± 0.051, 

= 1.761±0.013±0.070 and r^^) = 1.98±0.02±0.13, and at 80 GeV r = 1.548±0.008±0.041, 
r(^) = 1.834 ±0.016 ±0.088 and r^^) = 2.04 ±0.02 ±0.14, where the first uncertainty is statistical 



and the second systematic. These results are based on QCD analytic parametrizations [|I4|, ^ 
of the scale dependence of particle multiplicity in jets. We also present results which utilize 
polynomial parametrizations, or which assume no functional dependence for the growth of gluon 
jet multiplicity with scale, and obtain consistent results albeit with larger uncertainties. Our 
results are in agreement with the QCD prediction that r'-^-' should be closer to its asymptotic 
value of 2.25 than r^^\ and r^^^ closer to 2.25 than r, for the finite energies of our experiment. 
Our result for r^"^^ is the first experimental determination of that quantity. 



The results for r, r^^^ and r^^^ are compared to recent QCD predictions based on analytic 
and numerical [|13| techniques, as well as to predictions derived from the formalism of [jl5|. We 



find overall agreement between the experimental and theoretical results. 

Finally, we use the results on the energy dependence of Ng and Ng to determine an effective 
value of the ratio of QCD color factors, Cp^/C-p- We find Cx/C-p = 2.23 ±0.14 (total), consistent 
with the QCD value of 2.25. 
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